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We study the quantization of the electromagnetic sector of the Myers-Pospelov model coupled 
to standard fermions. Our main objective, based upon experimental and observational evidence, 
is to construct an effective theory which is a genuine perturbation of QED, such that setting zero 
the Lorentz invariance violation parameters will reproduce it. To this end we provide a physically 
motivated prescription, based on the effective character of the model, regarding the way in which 
the model should be constructed and how the QED limit should be approached. This amounts to 
the introduction of an additional coarse-graining physical energy scale M, under which we can trust 
the effective field theory formulation. The prescription is successfully tested in the calculation of 
the Lorentz invariance violating contributions arising from the electron self-energy. Such radiative 
corrections turn out to be properly scaled by very small factors for any reasonable values of the 
parameters and no fine-tuning problems are found. Microcausality violations are highly suppressed 
and occur only in a space-like region extremely close to the light-cone. The stability of the model 
is guaranteed by restricting to concordant frames satisfying f — |v maa; | > 6.5 x 10 _ . 

PACS numbers: 12.20.-m, 11.30.Cp, 04.60.Cf, 11.30.Qc 



I. INTRODUCTION 

The Myers-Pospelov (MP) model (l[ is an effective field theory that incorporates scalars, fermions and photons 
in a particle (active) Lorentz invariance violating (LIV) theory. It includes dimension five operators, together with 
the presence of a fixed time-like direction selecting a preferred frame. Such direction is assumed to arise from 
a spontaneous Lorentz symmetry breaking in an underlying theory and endows the model with covariance under 
observer (passive) Lorentz transformations. The modified free Lagrangian density is 

Cmp = ~4>* (d 2 + m 2 ) <P + (n^d^f - -F^F^ + X- (n"F^) (n a d a ) (n p e^ KA F KA ) 

+ (0„ - m) * + itf (n al a ) ( m + 77275) Kc^) 2 *, (1) 

to which we add the electromagnetic interaction via the standard minimal coupling. Such an effective theory is 
interpreted here as a model to describe the imprints at Standard Model energies of active LIV, codified by the 
dimensionless parameters p, £,?7i and 772, which is produced by drastic modifications of the space-time structure at 
a fundamental scale M, as suggested by some phenomenological models inspired upon developing quantum gravity 
theories 0, H, 0| and string theories @. Nevertheless, up to now there is no systematic derivation of a semiclassical 
approximation starting from a fundamental quantum gravity theory, for example, that could determine the exact 
nature of the possible corrections arising from such space granularity. This situation has prompted the construction 
and analysis of effective field theories models which capture the basic ingredients that we expect to survive at Standard 
Model energies. 

The additional Lorentz violating terms in JT]) are unique according to the following criteria: (i) quadratic in the 
same field, (ii) one more derivative than the corresponding kinetic term, (iii) being gauge invariant, (iv) being Lorentz 
invariant, except for the appearance of n^, (v) not reducible to lower dimension by the equations of motion and (vi) 
not reducible to a total derivative The model has recently been generalized to the non-abelian case including 
interactions arising from the fields associated to the Standard Model M. As such, it could be considered as a 
dimension- five-operator generalization of the Standard Model Extension 7] . In this work we will concentrate upon 
the simpler version of Ref . [l[ , particularly upon the proposed modified electrodynamics in its quantized version. The 
corresponding classical model has been thoroughly studied in relation to synchrotron radiation in Refs. |8[. Also, the 
self energy corrections of the model have been recently analyzed in [9( . Radiative corrections to LIV theories have 
been studied in Refs. [l(| and fine tuning problems have been discussed in Refs. [9I. flTj] . 

The point of view adopted in this work is to consider the quantum effective MP model (JXJ) plus the electromagnetic 
interaction as a perturbation of the Lorentz invariant theory, in the precise sense that after making zero the LIV pa- 
rameters encoding the corrections we must recover standard QED. Moreover, since all experimental and observational 
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evidence point to negligible LIV [12J , the radiative corrections arising from LIV should be accordingly very small. As 
we will sec in the sequel, this basic idea provides a guideline in the way one gives a meaning to the model, particularly 
in regard to its quantization and to the limiting procedure necessary to recover QED. 

Generally speaking, the dimension five operators make the theory of the higher order time-derivative (HOTD) type. 
This fact shows up in the Lagrangian ([1]) by the presence of third order time derivatives for the scalars, second order 
time derivatives for the fermions and third order time derivatives for the photons. It is well known that HOTD 
theories pose many difficulties for their implementation [l3j . [l4j], the most representative ones being the increase in 
the number of degrees of freedom with respect to the standard ones, together with the appearance of Hamiltonians 
which are not positive definite being unbounded from below. In this way, if one requires to treat the additional 
HOTD terms as a perturbation, a careful strategy is required. Fortunately, a systematic approach to carry out this 
task already exists in the literature [l5[ , [HI • 

In view of the above considerations a general strategy to define the quantum field theory extension of the MP 
model would be the following: (i) as usual, the starting point is the classical version of it given in Ref. [J. (ii) next, 
the application of the procedure in Ref. [IB] to the classical HOTD MP model would reduce it to a modified effective 
theory of the same time derivative character as classical electrodynamics. The procedure leads to field redefinitions 
plus additional contributions to the interactions, (iii) finally, this resulting classical theory would be considered as 
the correct starting point for quantization, which would be carried along the standard lines. The resulting quantum 
theory would then provide the basis for the calculation of interacting processes usin g th e perturbative scheme of 
quantum field theory (QFT). Some of these steps have been already carried out in Refs.fl71]. for the case of the scalar 
and fcrmion fields. 

Perhaps we should emphasize at this stage that we are dealing with two different classes of perturbations: the first 
one concerns only the LIV parameters, occurs at the classical level and serves to define the correct starting point 
for quantization. Once the resulting theory is quantized, the usual QFT interacting processes can be calculated, 
corresponding to the second class of perturbations. Both approximations should be made consistent when predicting 
a result to a given order in any of the LIV parameters. In this sense it is clear that we are not producing a quantum 
version of the full MP model, but only one which is adapted to our basic requirement of describing the LIV corrections 
as perturbations to QED. 

Since the model respects observer (passive) Lorentz transformations we consider the parameters p, £, 771 , 772 and 
M, to be invariant under them. Nevertheless, the general form of the four-vector describing the preferred frame 
is n M = 7(1, v), with 1/7 = vl — v 2 , so that highly boosted systems will greatly amplify the values of the LIV 
parameters which are strongly constrained in earth-based reference frames. Thus we also restrict the observer Lorentz 
transformations to concordant reference frames which move non-relativistically with respect to earth [l8j]. In the 
sequel we will give a quantitative characterization of such allowed observers. A further simplification is introduced by 
taking into account that the parameters p, ^, 771 and 772 are independent. In this way we set the field <f> together with 
the parameters p, 771 and 772 equal to zero. Then we deal with a minimal LIV extension of standard QED. 

The paper is organized as follows. In section II we discuss the classical MP modifications to electrodynamics. There 
we construct the corresponding Hamiltonian formulation in terms of canonically transformed fields that guarantee the 
appropriate normalization of the momentum squared terms in the Hamiltonian density, that includes the identification 
of the interacting sector. Section III deals with the quantization of the model in terms of standard creation-annihilation 
operators. The modified dispersion relations are identified and the Hamiltonian is shown to be positive definite for 
momenta k such that |k| < M/(2|£|). In Section IV we construct the modified photon propagator in the Coulomb 
gauge which is subsequently written in four dimensional notation by incorporating the static Coulomb contribution 
appearing in the Hamiltonian. Section V contains the physical motivation and specific proposal for our prescription 
that allows to understand the quantum MP model as a tiny perturbation of QED, according to the experimental 
and observational evidence of highly suppressed LIV. A coarse graining mass scale M << M is further introduced 
in the problem, dictated by the effective character of the model, and signaling the onset of the modifications in the 
space-time structure. In Section VI we set up the general structure the electron self-energy calculation including only 
the modified photon propagator (£ 7^ 0) interacting with standard fermions (771 = 772 = 0). The scale M is taken 
into account via a factor of the Pauli-Villars type, designed to act as the appropriate regulator in the QED limit. 
Also we perform a power expansion of the self-energy in terms of the external momentum and identify those terms to 
be subjected to scrutiny regarding their suppressed character and good QED limit in the next Section. The general 
strategy for their calculation is presented in Section VII and all the LIV contributions to order £ 2 are accordingly 
obtained. One of such calculations is presented in full detail, while we only write the results for the remaining ones. In 
Section VIII we present a preliminary study of the microcausality violation in the model by identifying the space-like 
region where it occurs, together with an estimation of the magnitude of such violation. The final Section IX contains a 
summary of the work. The notation and conventions are stated in the Appendix [Al which . together with the Appendix 
IBl contain information relevant for the specific calculations in the paper. In Appendix [C] the relationship between 
the modified photon propagator in different gauges is stablished. The last Appendix IDl includes the definitions of the 
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LIV contributions which calculation is not fully developed in the text. 



II. THE MODEL 



With the simplifications stated above we consider the modified the photon sector 

C-t = + X KF F ) (n a d a ) (n p e?» KX F KX ) - J»A U , (2) 

where the electromagnetic current J^ 1 will be subsequently realized in terms of unmodified spin 1/2 fermions, according 
to the choice 771 = 772 =0. Our general strategy will be first to quantize the photons and subsequently to consider the 
interaction, via the standard minimal coupling, with the unmodified quantum fermions. 
The equations of motion in the Lorentz gauge are 



(> A d 2 - 2g (n • d f n p eP^ x d a ) A x = J v . (3) 

In order to get a better control of the LIV modifications we find it convenient to work in the Hamiltonian scheme, 
so that we switch to a 3 + 1 canonical formulation of the problem. Taking advantage of the remaining observer Lorentz 
invariance of the model, we choose to work in the rest frame = (1, 0), where the free modified photon contribution 
is 

£7 = \ (A 1 + dtA ) 2 - I /•,,/•- + ge^ k A*d 3 A k - J»A„ g=t. (4) 

This choice has the advantage that, up to a total derivative, the resulting system is not of the HOTD type. Never- 
theless, it exhibits in a simpler setting most of the questions associated to the quantization of the full MP model. In 
addition, let us emphasize that we will carry the quantization without any approximation in the parameter g. 
The canonical approach gives the following momenta 

n = |^ = 0, U i =^=A i +d i A°+2ge^ k d j A k , (5) 
oA a 0A 1 

together with their Poisson brackets, 

{A%x),Il j (t,y)}=6}6 3 (x-y). (6) 

The next step is to construct the Hamiltonian density He = HiA* — C, which is 

He = TUA* - X - (A* + d z A°y + ^FijFV - gr* 1 '. Vi)J ! + J»A, L . (7) 

In order to write the velocities in terms of the momenta it is convenient to consider the combination Hi — diA together 
with the operator 

M ik = ^ik + 2ge l3k d 3 ) , (8) 

in the second Eq.Q. To solve for the velocities we need the inverse of the operator M %k for which we obtain the 
exact non-local expression 

(M-y = (1 + 4 1 g2y2) - 2ge^d r + Ag^d,) . (9) 

In this way we solve 

A* = (M- 1 )** (Tlj - djA°) , (10) 
which we substitute in Eq. ([7]). The result is 

He = (n p - d p A°) l - (M- 1 j pr (IL. - d r A°) - I (da ) 2 + \f io F^ + J»A,. (11) 
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Integrating by parts and using some of the properties for (M 1 ) written in the Appendix A, we arrive at 

He = \n p (M- 1 ) pr U r + (d p U p + J°) A + \f %3 F^ - J l A\ (12) 

It can be verified that the corresponding Hamilton equation of motion reproduces the correct expression p0[) for A 1 . 
The canonical variables can be written in the convenient form 

nf = A T + 2ge» k 8 j A T , Uf=A l L + diA°, A\ = ^di(d k A k ), (13) 

where we are using the standard definition for a transverse (T) and longitudinal (L) decomposition of a vector field 
U = + Ul, where V ■ Ur = 0, V x = 0. In the case of the velocities the separation leads to 



1 frii n^iria\nT ki 1 tt-L , a 1^9 



,2 



with the notation W = yl + 4g 2 V 2 . As in the usual case, Ao is a Lagrange multiplier leading to the Gauss law as a 
secondary constraint 

6 = dilk + J° = 0, (15) 

which can also be understood as arising from the time derivative LTo of the primary constraint Ho ~ 0. The evolution 
w leads to current conservation in such a way that we have only two first class constraints as in the standard 
case. In terms of transverse and longitudinal variables the Gauss law is written as 

d.Uf + J° = 0. (16) 

The equation of motion 

d l F l0 = +J°, 

yields 

A ° = -^2 (jo + dodiA*). (17) 
At this stage we select the Coulomb gauge by choosing 

A° = -^J°, Ho = 0, d k A k =0 - Ai(x) = 0, I$ = ~d i J . (18) 

The dynamical variables are contained only in the transverse modes 

Tlf = A t T + 2ge^ k d J A k , AUx), (19) 

which satisfy the Dirac brackets 

{AUt,x), Ill(t,y)}= U^-^U^-y). (20) 
Using repeated integration by parts in the Hamiltonian, together with the transversality condition, we arrive at 



Hc = J d* x QlLj (M-T ^ + \ J° (-^) J° + \F V FU JM t 



(21) 



Our final goal is to express the dynamical fields in terms of creation- annihilation operators, corresponding to mod- 
ified frequency modes, satisfying standard bosonic commutation relations that will reproduce the field commutation 
relations arising from the correspondence principle applied to the respective Dirac brackets. To this end it is necessary 
that the relation II = A holds, which is equivalent to require that the kinetic term of the Hamiltonian density be 
normalized as ^n 2 . In order to achieve this we perform the canonical transformation (At — ► At, II t — ► fl T ) given 
by 
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The non-zero transverse Dirac brackets for the variables A l T and IlJ are the same as in Eq. (l20|) in virtue of the 
canonical character of the transformation. Rewriting the Hamiltonian (|2ip in terms of the new variables leads to 

He = J dPx (jfijnj + \a\ [<T* - 2ge™*>d n ] A% + \j° (-^j J° - J*A^(A T )j . (23) 

Let us emphasize that in the last interaction term A % T is a functional of the dynamical field A J T . In this sense the 
electromagnetic vertex will be modified with respect to the latter field but will retain the usual structure with respect 
to the former. In this way, some care is required when implementing the perturbation theory starting from the zeroth 
order Hamiltonian written in terms of A l T and flj. 



III. THE QUANTUM THEORY 



Now we have the basic ingredients to proceed with the quantization of the modified photon field. We start form 
the usual plane wave expansion of the operator A l T (x) 



A^x) 



d 3 k 



£ 



1 



a A (k) e l (A, k)e- ik( - x > x + a[{k) e"(A, k)e + ' ifc ^' 



in terms of creation-annihilation operators a x (k), a A (k), respectively. The notation is 

[k(X)}^ = (w A (k), -k), fc(A) • x = u x (k)x° - k • x, 



(24) 



(25) 



where the modified normal frequencies will be consistently determined. The properties of the polarization vectors 
£*(A, k), A = ±1, chosen in the circularly polarized (hclicity) basis, arc collected in the Appendix B. The momenta 
are given by 



nfOr) 



d 3 k 



1 



2^ A (k) 



i-iux) a x (k) e l (A, k)e~ lk ^ x + (tw A ) a{(k) e"(A, k)e+ ife ( A )'* 



Assuming the standard creation-annihilation commutation rules 



a A (k), a{,(k') =S xx ,5 3 (k-k') 



and starting from (124[) and (|26[) we recover the basic field commutator at equal times 



[^(t,x),flT(t,y)] 



5 %3 



V 2 



S 6 (x-y), 



which is the expected result after the canonical transformation. The corresponding equations of motion are 

V 2 



K = d A 9 T , 



til 



Going to the momentum space we can obtain the modified dispersion relations from 

u; 2 A T (k) = |k| [A T (k) - 2. 9 zk x A T (fc)] , 
l-4g 2 \k\ 



which reduces to the diagonalization 



zkx — > Alkl 



when the vector potential is expressed in the helicity basis. In this way 

/ li |2 \ 



1- V|k|- 



[l-2A g |k|]A A , 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 
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yielding the modified energy-momentum relation 



(k) 



[l + 2A<?|k|]' 



(33) 



which is exact in g. With no loss of generality we assume from now on that g > 0. Let us notice that the four-vector 
[fc(A = +1)L is spacelike, while [fc(A = — 1)] M is timelike. At this stage we are confronted with two problems that arise 
rather frequently in LIV theories: (i) on one hand, the frequency w_(k) will become imaginary when |k| > l/(2g) 
and diverges when |k| = |k| maa; = l/(2g). From an intuitive point of view we consider l/(2g) as the analogous of the 
value |k| max = oo in the standard case and we will cut all momentum integrals at this value. The introduction of the 
coarse graining scale M << M, explained in more detail in Section V, effectively produces the more stringent and 
smooth cut-off 



5|k| <gM « 1. 



(34) 



(ii) on the other hand, since [k(X = +1)]^ is spacelike, we can always perform an observer Lorentz transformation 
such that w + (k) becomes negative thus introducing stability problems in the model. For a given momentum k this 
occurs for l/^/l + 2g\k\ < |v| < 1. Then, the condition (|34|) leads to the requirement that the allowed concordant 
frames in which the quantization will remain consistent are such that 7 < l/y / 2gM, with respect to the rest frame. 

Our next step is to verify that the resulting free (J M = 0) Hamiltonian is in fact positive definite and has the 
expected expression in terms of the previously introduced creation-annihilation operators. Let us begin with the 
kinetic term 



H 



KE 



d A x uf n 



Tfvr 



which leads to 



H 



KE 



WA(k) 



a A (k)a A (-k)e-^( k ) t 



— — a A (k) a A (k) 



h.c. 



(35) 



(36) 



in terms of the creation-annihilation operators. 
The potential term contribution is 



Hpot 

which analogously reduces to 

H POT = \ [d 3 k 



w 2 



[S rp - 2ge rnp d n ] A\ 



M(k) 



T- 



[a\(k) a\(-k) ] e 



-i2cjx(k)t 



^ a i(k)a A (k) + ^. C . 



(37) 



(38) 



Here we have made use of the dispersion relations 
final expression 



together with Eqs. (|B3p . (|B4|) . This leads to the expected 



H=\f rf 3 k^[a A (k)4(k) 



4(k)a A (k) WA (k), 



(39) 



arising from the cancelation of the time dependent terms and including the modified frequencies (|33[) . Thus the 
Hamiltonian is hermitian as far as the frequencies remain real, which is the case in the region |k| < l/(2g). 



IV. THE PHOTON PROPAGATOR 



In this section we calculate the free modified photon propagator starting from the definition 

l A lJ (x,y) = (0\T(Af(x)Aj(y))\0) 

= 9 (x° - y°) (0| Af (x) Aj (y) |0) + 9 (y° - x°) <0| Aj (y) Af (x) |0) 



(40) 
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where Ay (x, y) = Ay (a; — y) as can be seen from the expression 



(0|lf (x)IJ(2/)|0) 
Here we introduce the notation 



d 3 k 



T — 



-i(x-y)- k(\) 



(2tt) 3 Y 2w A (k) 



e<(A,k) 4(A,k). 



Fy(A,k)= ei (A,k) £*(A,k), 
which leads to the second vacuum expectation value in (|40p 

{0\Aj(y)AJ( X )\0) = J -^E^^-^F^U). 
We are interested in expressing the propagator 

Ay (*) = 



(2tt) 



4 e- ifc -Ay(fc), 



with z^ 1 = — yf 1 , in momentum space. To this end we start from the expression 
d 3 k 



iA i3 (z) 



Mk-z 



(2tt) 3 



E (A, £) + { - ZO ) g -^F* (A, -k) 



and introduce the standard representation 



(z a ) = i lim 



da 1 



e^o 7-oo 27r a + ie ' 
in order to calculate the corresponding Fourier transform. The result is 



A (M-W £i(A,k) £ ;(A,k) £j .(A,-k) e ?(A,-k) ' 
y V V 2w a (*o - + te) 2w A (fc + - «e) 



Using Eqs. (|B5j) and (|B6|) we rewrite the propagator in the form 



Ay (fe) = 



(*§-(• 



^A - «£) 



+ iX 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 



After performing the summations according to (|B8[) , (|B9[) and (|B10p we arrive at the following expression for the 
modified photon propagator in the Coulomb gauge 



Ay (fc) = 



1 



(fc 2 -4 9 



((fc2) 2 ~4 5 2| k | 2 fc4) 

Let us verify the correct limits when g = 0, where lo + 



2| k| 2 ^ 



(49) 



w~ = a; = |k|. In this case the first sum in the RHS of 
Eq. (|48[) gives the standard transverse propagator, while the second sum cancels out. 

We would like now to extend the above propagator, which is defined in the transverse sector, to the whole four 
dimensional space in such a way that the current-current interaction is described by 



d 4 k J tt (-k)A lu/ (k)J v {k). 



(50) 



This is achieved by incorporating in Eq. (|48p the Coulomb term appearing in (|2"3"|) in a manner analogous to that 
described in Ref. fijl. The final result is 



It, 



2^ {kl-LuS+ie) 



1 <5u(Ao ~ iXn p t 



p i»> i k i 



(51) 
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where we have reinserted the vector n p — (1,0). 

The last step in the construction is to perform the sums over A in (|5ip using the corresponding expressions in the 
Appendix B. The result is 



'//-< 



v (k 2 - 4 ff 2 |k| 2 fc 2 ) + Ag 2 \k\ 2 kl5^8 v0 + 2gn p e li r (ik a ) |kf 



([k 2 f- V| k | 2 fc4) 

The propagator obtained directly from the equations of motion ([3]) in the Lorentz gauge is 

1 



(52) 



A^(k) 



((k 2 



V| k | 2 fc4) 



'k 2 T\ illV + 2igk 2 e mr k m r]i^r] rl , 



k 2 



'-hkM + 



k 2 



-VOnVOu 



(53) 



In Appendix [C] we have calculated the propagator A^- corresponding to the fields A Z T starting from A^ given by 
and performing the canonical transformation (|22p . Moreover, the subsequent inclusion of the Coulomb term in Ajj 
leads exactly to the four dimensional propagator A^„ in (|53p . It is important to emphasize that the Hamiltonian 
(|2"3")l has a noninteracting sector described by the fields Aj,,f[J but induces an interaction density given by JiA l T — > 
iV^-A'*], where A 11 propagates according to ([5"3"|l. 



V. THE PRESCRIPTION DEFINING THE EFFECTIVE QUANTUM MODEL 

The main goal of this work is to study the possibility of defining the MP model as a perturbative extension of 
standard QED, that is to say as a model which continuously interpolates between a LIV theory and a Lorentz 
preserving one. This is to a large extent motivated by the very stringent experimental and observational limits set 
upon the parameters that codify such LIV. A construction exhibiting this interpolating characteristic has been already 
presented in Ref. (20| . but there the LIV was codified by a dimensionless parameter, as opposed to the situation here. 
As we will explain the sequel, the effective character of the model requires the introduction of an additional mass 
scale M that provides the analogous dimensionless parameter (gM). 

Another point that requires attention is the upper limit |k| max — l/(2g) set by the modified dispersion relations 
(|33[) . which guarantees the absence of imaginary frequencies together with that of a non-hermitian Hamiltonian. We 
consider these facts as indications of the effective character of the model. Assuming for a moment that 1/g ~ Eqq ~ 
Mpi anc k the above upper limit would mean that one is probing distances of the order of the Planck length, where we 
expect quantum gravity effects to be so important that the continuum properties of space might be no longer valid, 
thus invalidating the use of an standard effective field theory. This means that we need to introduce an additional 
coarse-graining scale M under which we can safely consider space as a continuum and apply effective field theory 
methods. Thus we require 

M «-, (54) 
9 

In this way, the upper limit |k| max = l/(2g) can be considered as a mathematical limitation in our model, analogous 
to |k| max = oo in the standard case. The physical limitation of the model is settled by the scale M and requires to be 
imposed by an adequate smooth regularization procedure that cuts down the corresponding degrees of freedom over 
this scale, which occurs a long way before energies of the order sa 1/g are reached. In this manner the relation (|54|) 
imposes a definite prescription to recover standard QED: (i) first set g — > for fixed M and (ii) then set M — ► oo. Let 
us emphasize that at the level of the effective model, the theory is finite and certainly will have an explicit dependence 
upon the physical parameters g and M. Now comes the question on how do we introduce the scale M. Intuitively we 
think of M as the parameter that will regularize the divergent integrals that will appear in the limit g — > describing 
standard QED. This suggests that we introduce this parameter via a Lorentz covariant smooth function X(fc), of the 
Pauli-Villars type for example, with the same characteristics that one would require in order to regulate standard 
QED. A natural choice for T(k) in our calculation of the electron self-energy is 

1 1 1 1 / M 2 \ 

k 2 m2 + ,/ ( fc ) = k 2 m 2 + ie - fc2 _ M2 + lt = fc2 _ m2 + ie ( M2 _ k 2 _ l€ ) ■ M»m. (55) 



In this way we are also imposing no additional LIV besides that arising from the original modifications to the dynamics 
encoded in the parameter g. 
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VI. THE ELECTRON SELF-ENERGY 



As a first step in testing the proposed construction we consider the calculation of the electron self-energy with 
the dynamical modifications introduced only via the LIV photon propagator. Let us recall that the perturbation 
theory based upon the Hamiltonian P5|) indicates that the photon propagates with A^ u given by ([55)1 . Moreover, we 
will focus upon the LIV contributions that could produce fine-tuning problems associated to the would be divergent 
contributions arising in the limit g — > 0. 



The starting point is 



d 4 k 



where we have introduced the scale M via 



(7 (p — k) + m) 



{{p — k) 2 — m 2 + ie) 



-."^.(AUIA-I^l-Ikl). (•)<!) 



M 2 



The ^-function is there to guarantee the reality of the frequencies a;.\(|k|) entering the calculation of the photon 
propagator in Section IV. Let us observe that the expression (|56p is finite. 

Next we find it convenient to expand the self energy in powers of the external momentum 



^> - ^ + Kw)^/ + 2 iPfJ,/'" +CMi (58) 

where each coefficient in the expansion is a matrix written in terms of some elements of the basis in the 4x4 
space of the Dirac matrices. We have considered up to second derivatives in the external momentum because the 
additional corrections to the numerator of the photon propagator (|53|) of order gk and (gk) 2 make those derivatives 
power counting divergent, as opposed to the QED case. The fact that we are violating Lorentz transformations in 
the boost sector, while maintaining rotational invariance would naturally split the above expansion into a time plus 
space structure. The expansion of the above coefficients in the gamma matrix basis will be denoted by 

-* _ u/src ( _ u/9 r c ( d£ g \ _ ma rC 



where we use the standard basis 

r c :T 4 = i, r" = y, r"" = o*", r 5 = 7 5 = ^VtV, r 5 -" = 7 V- (60) 

This allows us to rewrite the self-energy as 

£»(p) = (wg+jfWfao + \v> v W\^ T c + 0( P 3 ). (61) 

In order to deal with the calculation of such coefficients it is convenient to separate the modified photon propagator 
([55)1 into its even and odd parts 



A^(k) = A™ n (k) + A°f(k), (62) 
and rewrite them in the more compact form 

A^T(fc) = V^Fi + ^(|k| 2 W , ~ hkrS^l) = A— (k) = A™ e ™(-fc), (63) 
(k) = iF 3 k r e lmr Wm = _ A odd ( _ fc) = _ A °^(fc), (64) 

where 

F ^ F ^p/k 2 2gk 2 

({k 2 ) 2 -Ag 2 \k\ 2 ki) [(k 2 f-Ag 2 \k\ 2 kt) ((k 2 ) 2 - 4g 2 \k\ 2 kt 



are even functions of k and kg ■ 
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From the general expressions for the contributions in (|58[) , together with the symmetry properties of the propagator 
plus the symmetrical integration over the three-momenta it is possible to determine that the non-zero contributions 
to T, 9 (p) are 

= AI + iA 7 *y 7 fc e «* + p°B 1 ° - p i (CY - iC 7*7*6"*) + ^ (p ) 2 (oi + i/tyV 7V 3 '*) 

+ ip 2 [EI + iEYl j l k e ijk ) + iF p°p l (7°7 J 7 fc e yfc ) + 0(p 3 ) . (66) 

We will be interested in analyzing only those terms that could give rise to a finite and possibly unsuppresscd LIV 
contribution when g — > 0. In this limit we should recover QED, which is parity conserving so that we know that the 
electron self energy must have the form. 

£9=» = Wol + Wi ( P ^) + ^ P ^I + 0(p 3 ). (67) 

From this perspective all parity violating terms A, C, D and E in (|66p are subject to scrutiny and they should be 
finally suppressed. On the other hand, the parity conserving contributions can be rearranged in the following way 

Y? + (p) = AI+ (B - C)p° lQ + C(p" 7|t ) + ^±E1 (pO) 2 / _ + 0(p3)i (68) 

so that according to our prescription we expect 

Um(S-C)=0, \im(D + E)=Q, \im(A, C, -E ) = (W , W x , W 3 ). (69) 

g^O g^O g^O 

The general strategy to evaluate the required integrals is the following. The structure of the denominators D 
entering in them is of the form 

V= \(k 2 ) 2 - 4. 9 2 |k| 2 fc4 + ie] [k 2 - m 2 + ie] , (70) 

which can be rewritten 

V= (l - 4 3 2 |k| 2 ) [fc 2 - {oj 2 _ (k) -ie)] [fc 2 - K ( k ) -«)] [ko 2 - (E 2 (k) - ie)] . (71) 
Within the region of integration (|k| < 1/ (2p)), the poles in the complex fco plane have the form 

fcoi=£(|k|)-ie, fc 2 = -«?(|k|)+ie, (72) 



with £ (|k|) > 0. Here £ (|k|) stands for any of the involved energies u>± (k) and E (k) = Vk 2 + m 2 . In this way it 
is always possible to perform a Wick rotation to the Euclidean signature such that fco = *^4- Due to the remaining 
rotational symmetry, together with the symmetrical integration over k, one is finally left with only two integration 
variables which are and |k| that can be conveniently rewritten in polar form. 

VII. THE LIV CONTRIBUTIONS 

In this section we present a detailed calculation of the corrections W? y C to the electron self energy arising from 
the even sector of the photon propagator corresponding to the (B — C) term in Eq. (|68|) . The calculation of the 
remaining contributions goes along similar lines and we only give the final results. 

A. General structure of the contributions 

As a first step it is convenient to split them into the following temporal and spatial pieces 



<9£ \ 



7 o 2fc°(m — 7 fc) 



[fc 2 — m 2 + ie] [fc 2 — m 2 + ie] 



^A^(k)J(k), (73) 
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M =ie 2 



rf 4 fc 



Y 2k l {m~ 1 k) 
[fc 2 — m 2 + ie] [k 2 — m 2 + ie] 2 



^A^(k)J(k). 



(74) 



and further separate each contribution according to the even and odd pieces of the photon propagator A^ u (k). To 
simplify the notation we have introduced 



J(k)=l(k)> 



1 

2^ 



(75) 



1. The contributions from the even part of the propagator 
Taking the even part of the photon propagator in (|73p , the temporal component of the derivative is 



0EY 



iireven r M ,V,0„2 f j ^ 



1 



2&q 



(fc 2 — m 2 + ie) (fc 2 — to 2 + ie) 2 



where we have used 



^M^O^f A even o~,0 771 

7 7 7 A u„ = ~^7 fi- 



f-2Fi) ./(fc), (76) 



(77) 



The function i*i was introduced in Eq. (|65[) and from (|T6[) we see that the only contribution is given by the 
component WTqxq. Let us define the quantity 



1 



2/cq 



(fc 2 — m 2 + ie) (fc 2 — m 2 + ie) 2 



((fc 2 ) 2 - 4 5 2 |k| 2 fc 4 ) 



/(fc). 



The spatial contribution is (no sum over i) 



d P l ) P =o 



— ] rr/efen -pAf - 2 

7T- i - yy{i}M L 



-ie 



d*k 



1 



(fc 2 — m 2 + ie) (fc 2 — m 2 + ie) 2 



where we use 

/7"7"A;r = 2 7 * + (fc 2 - |k| 2 )F 2 ) 

The rotational invariance of the three-momentum integration leads to 



(fc 2 — m 2 + ie) (fc 2 — m? 



> 2 + ie) 2 \ (( fc 2) 2 _4 5 2 |k| 2 fc4) 



2|k| 2 /3 



/(fc). 



(78) 



(r^AZ en )J(k), (79) 



(80) 



(81) 



B. Calculation of the (B — C) contribution 



From Eqs. ([78]) and (@T]) we have 



u r : ^ I d ' k ^ + W fc2 " ^ |k|2fc " ^ - m2 + ^ /(fc) 

(2tt) 4 (fc 2 -m 2 + ze) 2 ((fc 2 ) 2 - V|k| 2 fc4) 



(82) 



In order to calculate the non covariant integrals of the above type, together with those in the Appendix D, we give 
some details of the procedure sketched at the end of the previous Section. Basically we implement the following steps, 
(i) First, we perform a Wick rotation to a Euclidean signature, such that 



fc =ifc 4 , fc 2 = -(fc| + k 2 ) = -k% <Tk E = i47r|k| z dk A d|k| 



(83) 
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(ii) Second, since we are maintaining rotational invariance we are left with only two variables 

1 



oo < k 4 < +00. < Ikl < 



2.9 ' 



In this two-dimensional space we introduce the following polar coordinates 

ki = rcosa , Ikl =rsina, 



(84) 



(85) 



where k 2 E — r 2 . Next we have to integrate over the rectangular strip defined by (f8~4")l and we choose first to integrate 
over r and subsequently over a. In this way we have 

d^k — i J d^kg = i 4tt / da siir a 
Applying the above procedure to Eq. ([82|l we have 



l/(2fl£ 



r A dr. 



(86) 



B — C = -Ae J 



d 4 k E 


^E i^OE ~ Ik 2 ) + 


I^IMXe (-fcl-m 2 + ||k| 2 )" 


M 2 


(2,f 


(k E + to 2 


f ((fc £ 2 ) 2 -V|k| 2 ^ £ ) 


(M 2 + k 2 E ) 



^-Ikl 



(87) 



Introducing the polar coordinates (|85|) yields 

2 /--it /■l/(2gsina) 



^ 3 ./o 



da sin 2 a 



(cos 2 a ~ \ sin 2 a) + |<? 2 sin 2 a cos 4 a{— r 2 — m 2 + |r 2 sin 2 a 



(r 2 + to 2 ) (l — 4g 2 r 2 sin a cos 4 a 



M 2 



(M 2 + r 2 )' 
The required radial integrals are 

K x {a) 



M 2 



l/(2 9 sina) r 3 

(r 2 +m 2 ) 2 (l-r 2 c 2 (a)) (M 2 + ) 



dr. 



(89) 



^ 2 (a) 



1/(2^ s i n Q ) 



M 2 



dr. 



which can be exactly calculated, yielding 



M 2 



(c 2 to 4 + M 2 ) /A 2 
p 2 (A 2 ) 2 n 



(r 2 + to 2 ) 2 (1 - r 2 c 2 (a)) (M 2 + r 2 ) 

c 2 ln(l - c 2 b 2 ) 



M 2 /A 2 



(A 2 ) 2 ^ V M2 



p 2 77 



pA 2 „A 2 



(90) 



(91) 



#2 



M 2 



(2to 2 M 2 + to 4 (c 2 A/ 2 - 1)) /A 2 
In 



p 2 (A 2 ) 2 



M 4 . M 2 ^ ln(l-c 2 ^ 
(A 2 ) 2 7y ^Af 2 



b 2 m 2 



P 2 V 



P A 2 m A 2 



The notation is 



A 



p = 1 + mV(a), 77 = 1 + MV(a), A = m — M , c(a) = 2gsinacos 2 a, 

2 , l2/ \ A 2 _ „. f 2 .2/ - 1 



m 2 + b 2 (a), A J M = M 2 + b 2 {a), b{a) 



2g sin a' 



(92) 



(93) 



In order to simplify the results by including only the dominant terms, we will expand the above expressions in powers 
of g 2 . This is justified since the expressions ([9Tj) and l[92|) are free of poles. Up to order g 2 , the remaining integrals 
over a will be of the form 



sin p acos 9 a, sin p a cos 9 aln(l — cos 4 a), sin p a cos 9 a ln(sin a) , 



(94) 



with p, q integers. These integrals contribute only with finite numerical factors, which are not very relevant in order to 
establish the correct QED limit of the LIV terms and only the final numerical results will be presented. Nevertheless, 
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we will isolate the exact g 2 independent contribution and we will show that the angular integration produces a zero 
contribution, thus eliminating any indication of fine-tuning. In all the remaining contributions proportional to g 2 we 
will further expand in powers of m/M and retain only the dominant terms. In this way we will need the approximate 
expressions 

/ jy4 f M\ M 2 \ ( \ 

^ = [ (m2 _ M2)2 ln U J + 2(m 2 W) ) + 2{gm)2 Sln2 a C ° S4 a [ 1 + A H2gm SlnQ) J 

-2(gM) 2 sin 2 a ^1 + cos 4 a (2 In (2gM sin a) - cos 4 a In (l - cos 4 a))^j , (95) 

g 2 K 2 = -(.9Af) 2 ^ln(2M 5 sina) + iln(l - cos 4 a)^j + ( 5 m) 2 ^2 ln(2. 9 m sina) + ^ . (96) 

It is important to observe that the exact g 2 independent term, contained in the first bracket of Eq. (|95|) gives a zero 
contribution in virtue of the angular integral factor 

sin 2 a (cos 2 a — — sin 2 a) da = 0. (97) 
o V 3 / 

Performing numerically the remaining angular integrations in the proposed approximation we obtain 

B — C = ^| (gM) 2 ( - 0.070 + 0.010 1%M)) 

-(gm) 2 (o.016 + 0.021 \n(gm) + 0.031 ln (^) ) | • ( 98 ) 
The remaining contributions from the even sector are 

(gM) 2 (0.75 + 0.047 ln(#M)) - (5m) 2 (o.014 + 0.047 h%m)) J , (99) 



e 2 m M 2 f\i 

~ ~ I 2(m 2 - M 2 ) 11 



^( ln S)-l)' « 



^ + S = -^' ' ' ' 



Finally, the odd contributions are 

e 2 f . ,0 / , . ,A o/ , ,\\ ^ e l , , I '. f m\ 1 



A= ^g M 2 (0.018 + 0.063 lnfoA/)) - m 2 (o.026 + 0.063 ln(.gm)) , C = ^(gm) hn(^) +-j , (101) 



5ge 2 / fM\ 4\ - ge 2 /13 f M\ 17, 

^^MmJ+s)' S ^(24 ln b) + 32 )- < 102 > 



The results obtained above, in the framework of our prescription to recover QED, have precisely the expected property 
that reduce to zero when we turn off the LIV correction parameterized by g, keeping M fixed. Also, the results are 
consistent with the fact that the unsuppressed contribution which we still expect to diverge even after we set g = 
and subsequently M — > 00, comes in the term A written in Eq. (|99[) . This term corresponds precisely to the mass 
renormalization contribution in standard QED. 



VIII. MICROCAUSALITY VIOLATION 



In this section we provide an estimation of the microcausality violation associated to our model. A comprehensive 
study of such violations is out of the scope of the present work. Microcausality violation has been previously studied 
in the fermionic sector of the Extended Standard Model, for example [l8| . 
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We work directly in the Coulomb gauge associated to our reference system where = (1,0). We only consider 
points x and x' which produce a space- like interval (x — x') 2 < 0. Unfortunately we cannot perform a passive Lorentz 
transformation to reach the system where xo — x' = 0, which might simplify the calculation. This is because such 
transformation will change n M = (1,0) into = 7(1, v) and then our system will turn out to be manifestly of the 
HOTD type, thus requiring the application of the perturbative process of Ref. (16J, which we have avoided in our 
particular reference frame. 

Even in the standard QED case there is a drawback when working in the Coulomb gauge, which is basically due 
to the apparent causality violation of the theory arising from the instantaneous character of the scalar potential. 
When dealing with the commutator [A^(x), A J T (x')], which is the naive starting point to test microcausality, this 



D(x — x 1 ) where D(x — x') is a function 



problem shows up because this commutator is proportional to ( — 

that has support only in the light cone. Nevertheless, the operator 1/V 2 , which is just a shorthand for the Green 
function l/\r — r'|, acting upon D(x — x') produces non-zero results outside the light-cone, thus yielding an apparent 
violation of microcausality. The canonical way of dealing with this problem is to calculate the commutators of the 
gauge invariant fields E and B for space-like separation. We will follow the same route here and we will discuss only 
the commutator 

[Uf (x), nj(x')] = -dl [DV(x - x')] = W(x - x'), (103) 
which is the analogue of the electric fields commutator in standard QED, with ILf (x) been gauge invariant. Here 

D ij (x - x') = [Air(x), A^(x')}. (104) 
A direct calculation starting from Eq.(32) leads to 

d 3 k 

(2^)3 ^ 2c^ A (k) 

and we use z M = (x M — x M ') in the sequel. Using the relations (|B5[) and (|B6[) from the Appendix we arrive at 

d 3 k 1 



d 1 Hx-x') = J ^ Ea^?k) ( e< ( A ' k y*( A ' k ) e ^ (x "^ ) - E< *( A . k ) e, '( A . k ) eifcx(a ~ x ' ) ) . 



(105) 



which can be rewritten as 
D l3 {z) 



E 



s ij 



(2tt) 3 2w A (k) 
d 3 k 1 



1 




iX 












ijm 




2 


|k| 2 _ 


~ y 


|k|. 





(2tt) 3 2^ A (k) 



1 




i\ 


ijm * 




2 


|k| 2 _ 


+ y 


|k|_ 





-i(u\z — k-z) 



i(w A z -k-z) 



(106) 



QiQj 



E 



1 



d 3 k 



d 3 k 

W 2|w A (k)| 
A 



(2tt)3 2|k|| WA (k)| 



-i(w\z — k-z) _ i(u>\z -k-z) 



-i(w A z -k-z) _ i(w A z -k-z) 



(107) 



Let us remark that this expression contains the correct limit when g = 0. In this case ui+ = oj+ = |k|, so that 
the contributions of each term in the J2\ are the same. After the summation, the first line of ()10T|1 reproduces the 
definition of the standard function D(x — x'), while the second line is proportional to J2\ A = 0. 
Starting from (|103|) yields 



n ij (z) 



d 3 k v/r+2A4k| r e -<( U> »-k..) 



g [W] 



(2tt) 3 

1 



E 



2|k| 

d 3 k 

(2^)3 



v/1 + 2A. 9 |k| 



O i(uxzo — k-z) 



-i(w A z -k-z) _ i(w A z -k-z) 



(108) 



where we have rearranged the above expression in such a way that the first line of (| 1 08[) recovers the standard QED 
result in the limit g — > 0, while the second line is equal to zero. In this way the microcausality violation is encoded 
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in the functions 



- ^E 



d 3 k v /l + 2A 5 |k| 



E 



(2tt) 3 

ri 3 k 
(2^)3 



2|k| 



V / T+2A#I 



e -i(^AZo-k-z) _ e i(o) A zo-k-z) 



g -i(wA2o-k-z) e i(w A zo-k-z) 



(109) 
(110) 



which are now acted by local operators only. 

Since we expect microcausality violations, we will estimate their impact arising only from the function V\. Notice 
that V\ (z) = —Vi(—z) as can be seen from the expression (| 109[) . After performing the angular integrations we obtain 



1 1 



2(2tt) 2 



-E 



V 2 9 



dk 



VI + 2Xgk 



[e lkr 



kr ikrl 



(in) 



where k = |k| and we have enforced the upper limit l/2g in order to have real frequencies w\(k) according to Eq. 
(I33p . The spacelike character of the interval is written as — r < zq < r. To proceed we introduce the phases 



$u(fc) = k[r 



1 



\/TT2Xgfc 



z , ®2\{k) =k[r + 



1 



VI + 2Xgk 



2() 



in terms of which we rewrite V\ as 

v, 1 11 



2(2tt) 2 ir 2 



x /-l/2fl 

/ dfc v/l + 2A 5 fc 



g'*lA _|_ g-**lA _ g**2A _ g- i$ 2 



(112) 



(113) 



In order to make an estimate of the region where microcausality violations occur we concentrate in the calculation of 
the momentum integrals appearing in Eq. (|113p . We apply the stationary phase method to the generic integral 



h = 



l/2fl 



dk f x (k) e -^ k \ f x (k) = a/1 + 2Xgk, 



where the relevant phases are given in Eq. Q112p . The general result for such integral is 



h = fx(k)e 



i* A (fe) 



1/2S 



dk 



(fc-fc) 2 



(114) 



(115) 



where k is the momenta that makes de phase stationary within the interval [0, l/2g]. 

We illustrate the calculation for the case of The remaining cases are completely similar and only the final 

results are written. The exact expression for the momentum k that extremizes $ia is given by the equation 



r 

Z() 



1 + Xgk 



(1 + 2Xgk) 2 

Observe that k appears always in the combination gk so that the solution will be of the form 



k = -x — , 

9 \ z o, 



(116) 



(117) 



where x solves the corresponding equation obtained from (|116|) . This is a complicated function of and to 
make some analytical progress the following approximation is made. We found that in the range of = 1 + e with 

e « 1, the exact curve x [J^j m we ll approximated by the straight line 



hx = - T - --1 , ->l, 

2g \z J Zq 

resulting from the expansion of the phase to order k 2 in Eq. (|112|) . which is 

<f> 1A (fc) = kr - (k — k 2 Xg) z Q . 



(118) 



(119) 
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This means that we are considering a space-like region close to the light-cone such that 

(1 - e)r < \z \ < r. 



(120) 



A posteriori we will verify that our results in fact fall within the range of the approximation. For this purpose it 
is convenient to rewrite the condition (|120|) by stating that the maximum allowed fractional deviation j^r 11 ! has to 
satisfy 



Az 



< e. 



(121) 



From now on it is convenient to separate the cases according to the sign of zq. For zq > the extremum (II 18[) has to 
satisfy the condition 



< fci 



A / r 



A 



1 < 



1 



2ff \z 

We observe that we have no solution for A = +1. The choice A = — 1 requires 

r 



In this way we have 



t 1 I r 

2ff \z 



< Z Q . 



^l-(fcl-) 



l + 2\gk 1 _ = J2- 



zo 



rf 2 $i^(fc) 
dk 2 



( r - z o) 

kgza 

= -tgzo- 



2 < z o < r , 



k—ki 



The case z < produces 



with 



<S> 1 x(k)=kr+(k-k 2 \g)\z \, 



2g \ \zq 



The condition 



0< k[> 



A / r 



2g Vko 



l < 



2.9 



(122) 
(123) 

(124) 
(125) 

(126) 
(127) 

(128) 



cannot be satisfied neither for A = — 1, nor for A = +1. The former leads to negative k' ix , while the latter requires 
t-^— r < 0. In other words there is no solution for z n < 0. 



The case of <S>2\(k) has solution only for z < and A = — 1. The results are 

x 2 



to- = 



1 

2^ Vl z o 



1 



'l + 2Agfc 2 _ = 
Substituting in (|113[) yields 

^(z) = 



$2-(fc2 

d 2 $ 2 _(fc) 



dfc 2 



= -2.9M 



< |z | <r, z < 



(129) 
(130) 



k=k 2 



1 1 

4(2tt) 2 ir 
1 1 



0(z o ; 



2z n - r 



Z{) 



' r 1/23 

g 4 9*0 



4(2tt) 2 


IT 


1 


1 


4(2tt) 2 


ir 


1 


1 



*(*<>; 



4(2tt) 2 



-zo, 



2z - r 

zo 

2z + r 
zo 

2z a + r 
zo 



r 1/23 



dk e+^Mk-k!-) 2 



dke +i9Z °^ k ~ k2 ^ 2 



2 ,l/2 fl 



dk e -*9Mk-k 2 -) 2 



(131) 
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where we can verify that V%(z) = —Vx(—z). Though this will not be relevant for our estimation of the microcausality 
violations, we can go one step further and estimate the remaining integrals in the following way. Introducing the 
change of variables u = Jg{k — we obtain 



/id 



V 2 9 - I 

dk e^aMk-k!-) 2 _ 

\/9 J-^/gki 



du e 



Substituting the value of results in 



(132) 



h± = 



V9 



2V9 \ 



2n/9 V z 



du e ±lz ° u 



1 

7a 



— I due ±lZ0U = 



2gz 



(l±z). 



The expression for 



h± = 



1/2S 



dk e ±i 9 z o(k-k 2 -) 2 _ 



1/23 



dk e ™M(fc-fc2-) 2 



can be obtained from (|133p changing zq by |zo|, so that we obtain 

h± = 

Then we have 

Vi(z) 



2g\z \ 



(IT*)- 



(133) 



(134) 



(135) 



1 



4(2ir) 2 ir 

1 1 
"4(2^)2^ 



2z 



2gz 



e 4 shi [1 — i) + e 4gz '> (1 + i) 



0(-zo) 



2z 



zq 



2g\z \ 



.-(■■+^o)^ 



:(l+£0) 



(1 - i) + e + 4 ^o (l + i) 



(136) 



Next we analyze the regions where microcausality is violated and provide an estimation of the amount of such 

±i (r-i» l) 3 

violation. In our approximation such violations occur when the functions e i3Z o do not oscillate rapidly enough 
to make V\ (z) equal zero in the space-like region. Thus we take the condition for having microcausality violations to 
be the region where the phases change slowly, that is to say where 



(r-\zo\Y 



< 1, 



(137) 



in which case the oscillations are very much suppressed. Let us concentrate now in the case zo > (the case zo < 
can be discussed in a similar way). We first examine the curves that limit the region of interest by considering the 
equality in Eq. (|137p . For a given r, the solutions of such equation are 



= r + 2g + 2^/g 2 +gr, z - = r + 2g - 2^g 2 + 



gr. 



(138) 



We observe that zq+ is always above the line zq — r, while zq_ is always below. Also notice that both curves tend to 
the light cone when g — * 0. The condition (1137|) is satisfied when 



2\J g 2 + gr -2g] < z < r, 



(139) 

because this region includes the case zo — > r which clearly satisfies (|137[) . That is to say, (|139[) determines the space- 
like region where Vi(z) is not zero, thus leading to microcausality violations. For a given r, the range of zq within 
that region is given by Azo = r — zo- ■ Then we can quantify the maximum time interval for which such violations 
occur by 



Az Az 1 / y-j— — 

= -[2y/g 2 +gr- 2g 



The expression in the RHS of (|140jl is a monotonically decreasing function of r with the following end points 



"Az " 



= 1, 



'Az " 



. 



(140) 



(141) 
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ro [cm] 


|Az |/r < 


1 


6. 


3 x 10" 22 


10 -io 


6. 


3 x 10" 17 


1Q -22 


6 


3 x 10" 11 


itr 27 


2. 


x 1(T 8 


1Q -33 


2 


x 1(T 5 



TABLE I: Upper bound on fractional microcausality violation [Azo| /f for distances r > ro 
That is to say, for the whole region r > r we can guarantee that 

r 

Recall that g = f/M, where £ is bounded by 10~ 10 when we choose A? = A/ P = 10 19 Gev (L P = 10~ 33 cm) [22| . 
Thus, taking g = 10 -43 cm and always considering the region r >> g, where we can trust the effective theory, we 
make some numerical estimations of the relation (I142p , which are given in Table [T] The third and fourth values of 
Table U correspond to distances given by ro = 10 11 Lp and ro = 10 6 Lp, which set a lower limit beyond which space 
becomes granular, according to the models considered in Refs. [ij and [23[ respectively. The calculated microcausality 
violations in Table I fall comfortably within the range determined by (|12ip required for the approximation to order 
k 2 in the phases (| 11 61 ) to be correct. 

IX. FINAL REMARKS 

In this work we have proposed a consistent quantization of the electromagnetic sector of the Myers-Pospelov (MP) 
model, [l| coupled to standard fcrmions, such that it can be realized as a perturbative correction of standard QED. 
By this we mean that in the limit where the Lorentz invariance violating (LIV) parameter g = £/M goes to zero one 
should recover the same quantum corrections arising in QED. Even though this sector of the MP model is not of the 
higher order time derivative type, up to a total derivative, some subtleties appear in the quantization of the photon 
field. The correct perturbative prescription is achieved by recognizing the effective character of the model via the 
introduction of a coarse graining scale M << 1/g, under which we assume that space retains the usual attributes 
which allow the construction of a standard effective field theory. Such cut-off scale is incorporated, in a smooth way, 
by means of a Lorentz covariant function of the Pauli-Villars type, which plays the role of a standard regulator in the 
QED limit and makes sure that all LIV is codified in the parameter g. The mathematical translation of this physical 
picture amounts to the following prescription in order to properly recover QED: first take g = 0, for constant M , and 
subsequently set M — > oo. The prescription has been tested in the calculation of LIV contributions arising from the 
electron self energy, which indeed provide the expected results. In this way the fine tuning problems found in Refs. 
0, EH disappear and one in fact recovers the correct zero limit for all the LIV corrections, which are indeed shown to 
be very small perturbations in accordance with the experimental and observational evidence. 

Some comments regarding the plausibility of the scale M in relation with the very stringent constraints already 
found for LIV are now in order. The combinations of parameters g = £/M, iji^/M, denoted collectively by 3/M, 
appearing in Eq. (TT]) are considered as remnants of a more fundamental quantum gravity (QG) theory, which include 
effects that make space no longer describable in terms of a continuum. Such parameters could arise, for example, 
in the process of calculating expectation values of well defined QG operators in semiclassical states that describe 
Minkowski space-time, which would be necessary to derive the exact nature of the induced corrections to standard 
particle dynamics at low energies. Let us emphasize that what is bounded by experiments or observations is the 
ratio S/M, so that a neat separation of the scale M and the correction coefficients 3, which could even be zero if no 
corrections arise, is not possible until a semiclassical calculation is correctly performed starting from a full quantum 
theory. Initially, the naive expectation was that taking M = Mpi anc k will be consistent with 3 values of order one, 
which is certainly not the case. Nevertheless, we should not rule out rather unexpected values of 3 or M until the 
correct calculation is done. 

Let us assume that we have identified the correct separation in 3qg / Mqg consistent with the experimental bounds 
and arising from a correct semiclassical limit of the QG theory. Then we will interpret Mqg as the scale in which 
quantum effects are manifest and where space is characterized by strong fluctuations forbidding its description as a 
continuum. Nevertheless, another scale M naturally should arise in this approach, which is the one that separates the 
continuum description of space from a foamy description related to quantum effects. That is to say, for probe energies 



Az 



(142) 
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E « M we are definitely within the standard continuum description of space where effective field theory (EFT) 
methods should apply. For probe energies E >> M we enter the realm of quantum gravity and there we assume that 
any EFT has to be replaced by an alternative description. It is natural that a very large number of the basic quantum 
cells of space characterized by the scale (1/Mqg) 3 will contribute to the much larger cells characterizing the onset of 
a continuum description, so that we expect M <C Mqq. 

The maximum allowed momenta |k ma;r | w Mqg /Sqg in the theory will be mathematically dictated by the positivity 
of the normal modes energies, Eq. (|33p in our case, and certainly constitutes an extrapolation of the EFT that can be 
considered as the analogous of taking the maximum momentum equal to infinity in the standard QED case. That is 
to say, we need to introduce an additional suppression of the excitation modes in our EFT which will be settled by the 
scale M, thus defining the effective energy range of the model. This is required by the EFT description of excitations 
in space which demands that the Compton wave length l/|k| of the allowed excitations be larger than the scale 1/M 
setting the onset of the continuum. The implementation of this proposal is directly related with our demand that the 
quantum model constructed from the MP theory be such that it produces a continuous interpolation between those 
physical results including 5^0 corrections and those predicted by standard QED (S = 0). In order to achieve this 
we have proposed the prescription fully described in Section V. 

Let us now discuss whether or not an estimate of the order of magnitude of the scale M in relation to Mqg makes 
sense. In our specific case the LIV contribution to the electron self energy produces an additional dimension four 
contribution to the Lagrangian given by 

e 2 - 

AL = 6^~f°id y, (143) 

arising from the (B — C) term in Eq. (|68p . Our calculation leads to a prediction dominated by 

\5\ ~1(T 2 x {gMf\\n{gM)\, (144) 

according to Eq. f|98[) . On the other hand, starting from the correction (|143[) together with bounds from the anisotropy 
of the inertial mass, the authors of Ref.0] have established the experimental bound 

\5\ < 1(T 21 . (145) 

In this way, we expect that the scale M is bounded in such a way that the theoretical correction (|144l) is much less 
than the experimental bound Q145p . that is to say when 

gM = f^L = 1CT 10 ^ << 0.65 x 1CT 10 , — > M«0.65M. (146) 

The above shows that it is safe and consistent with present observations to define a scale M much below the quantum 
gravity scale M. Proposals for additional scales M significantly smaller than M = Mp, that can be understood as 
signaling the transition between the standard space-time and that associated to the quantum gravity phase, already 
exists in the literature, [il. |23T|. 

Next we comment upon the behavior of our result for the electron self-energy under different momentum routings. 
For arbitrary internal momenta, the basic expression (|56p can be rewritten as 

E» = -ie 2 J ^^S F (k a + k^A^(-k a + k°)l(-k a + fcf) 9(l/(2g) - |k x - k|), (147) 

where k" — k a (k a + k% ) is the internal photon (electron) momentum respectively, with kf and k% being arbitrary 
momenta satisfying the conservation k% — p a — kf . Here Sp(k a ) denotes the standard fermion propagator. Since the 
integral (|147p is finite we are allowed to make the change of variables k a — fcf — > k a , which reduces the integral to 
the form (|56p and shows its invariance under momentum rerouting. 

The stability of the model is guaranteed by restricting the observer Lorentz covariance to concordant frames 
characterized by boosts factors up to 7 = \j\f2g~M. Using the bound (I146P the maximum allowed boost factor is 
Imax — 8.8 x 10 4 , which corresponds to a maximum relative velocity such that 1— |v maa; | > 6.5 x 1CP 11 . This condition 
certainly includes concordant frames that move non-relativistically with respect to earth. 

We have made a preliminary estimation of the microcausality violations in the model by looking at the commutator 
of two gauge invariant momentum operators (which are the extension of the electric field operators in standard QED) 
for space-like separation r > zq. The value of the corresponding function has been calculated using the stationary 
phase approximation and the condition for having microcausality violations requires that the exponentials oscillate 
very slowly. This means that the associated phases should be of order one or less, which defines a space-like region 
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extremely close to the light cone, rapidly approaching to it when the LIV parameter g — > 0. For a given value of r the 
width |Az | of such region is calculated. The fractional value (| Azo\/r) max which sets the upper limit for the allowed 
microcausality violation is subsequently estimated, leading to a typical value of \Azo\/r < 6.3 x 10~ 17 for distances r 
larger than the Compton wave length of the electron. 

In this paper we have studied the construction of the quantum MP effective model emphasizing the recovering of the 
correct QED limit in relation with the absence of fine-tuning problems. A summary of our results has been presented 
in Ref. [24j . Within the restrictions imposed we have established the basis of a sound perturbative scheme to proceed 
with the calculation of additional radiative processes. We defer for further work the analysis of the predictive power 
of the model in relation to LIV corrections to physical observables. 

APPENDIX A 

The conventions used in this work for the Dirac algebra are those of Ref. [2l| and we take K = c = 1. Also we have 
V, = diag(l, -1,-1,-1), 6 ij =+l,i=j, <5« 0. i t j. 

e 0ijk = e ii fcj e 0123 = £ 123 = = +lj 

A = (A^-Ai), k=(k i = -k i ), v = f^ 7 = a f V 



d r -ikr, dt^-iko, A • B = A l B\ (A X Bf = e ijk A j B k , (VxB)' = e'^/. (Al) 

In addition let us summarize some useful properties of the operator (Af — introduced in Section HT1 

(M-Y (M-y = 1 [(1 - 4.g 2 V 2 ) 6* - ige^d p + 4. g 2 (3 + 4. g 2 V 2 ) d 3 d r ] , (A2) 

(l + 4g z V z ) 

(Af-i) pl e ijkg. (M" 1 )^] = (1 - 4 5 2 V 2 ) e^dj + Ag (d pr V 2 - d r d p ) , (A3) 



(1 + 4<7 2 V 2 )' 

( M ~r = d + v V2) (r + v m-) - (A4) 

d j (M-y = (1 + 4 X g2v2) (1 + 4.g 2 V 2 ) ft = ft, (A5) 

(m-T nj = {l + l g2v2) (s- + 2 9 ^d r ) nj. (A6) 

APPENDIX B 

Here we present the properties of the polarization vectors in the hclicity basis (A = ±1), which are used in the 
expansion of the photon field in Eq. (|2"4")) . They satisfy the identities 

e*(A,k) = e(-A,k), k-e(A,k)=0, k x e(A,k) = -iAe(A,k), (Bl) 

£ *(A,k)-e(A',k)=5AA', e(A,-k) = e(-A, k), e(V, -k) • (k x e(A,£)) = -iA|k|5 AA /, (B2) 

e r (X,k)[6 r P + i2 9 e rm Pk m }eP(X',-k) = [1 + 2<?A|k|] S xy , (B3) 

e r *(\,k)[6 rp -i2ge rmp k m ]e p (\',k) = [I + 2g\\k\]5 xy . (B4) 
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The following combinations are useful in the construction of the corresponding propagator 



£ *(A,ky*(A,k) = \ 


lkl 2 _ 


-4 


|k|. 


e l (A,-kV*(A,-k) = 1 


{ Ikl 


+4 


ijm * 

|k|. 



together with the sums 



^2 £* ( A, k)e J * ( A, k) = ^ e j (A, — k)e" ( A, — k) = S ij 



±A ±A II 

The final calculation of the propagator in Eq. (|5ip requires the calculation of the following sums 

2(k 2 - V|k| 2 fcf 



E 

A 

E 



i 



A (^o-^A + i£ ) " (k 2 f- V|k| 2 fc4 + ze' 

2k 2 fc 2 



A 



A 



E 



(/eg - + ie ) (P) 2 - V|k| 2 fc4 + ie ' 
A 4 5 |k| 3 



( fe Q - W A + "0 



( fc 2) 2 „ 4 5 2| k |2 fc 4 + z£ - 



(B5) 
(B6) 

(B7) 



(B8) 
(B9) 
(BIO) 
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The purpose of this Appendix is to compare the propagator obtained directly from the equations of motion in 
the Lorentz gauge and further expressed in the Coulomb gauge, with the propagator (|40[) obtained directly in the 
Coulomb gauge after the canonical transformation (f22|) is made. 

From the equations of motion ([3]) in the Lorentz gauge we identify the momentum space operator 



0"*{k) = -fcV - 2gik 2 e 0l " T<l >k !J , 



(CI) 



which propagator is 



(fc2)2 _ Vfc 4 | k |2) 



2 2m Ag 2 k A 

~k T)iiv 4" 2ig/CQ£ P 9 kmTjp^Tjqp — kpkqS^S^, 



4<7% 4 |k| 2 



VonVov 



such that O^^A^p = 6". Separating the instantaneous Coulomb contribution 

J l (-/c)A„(fc)J J (fc) = J^(-fc)A^(fc)J iy (/fc) - J°(-fc)-L j0( fc ) 

|k| 

and using charge conservation [l^ | , we find the corresponding propagator in the Coulomb gauge 

1 



Ay(fc) 



ikr((fc 2 ) 2 -4^| k | 2 ) 



k ^|k| <5jj kikj 



2 5 fc 2 |k|V m ^ 



(C2) 



(C3) 



(C4) 



The Coulomb gauge propagator corresponding to the canonical transformation (|22|) and which was directly constructed 
from its vacuum expectation value definition (|40|) leading to the final result Eq. (I49|) is 



A,, = 



k| 2 (fc2) 2 - 4 5 2fc 4| k | ; 



(fc 2 - 4 5 2 |k| 2 /c 2 ) (|k| 2 % - kikj) + 2g\k\ 4 ie tmj k 



(C5) 
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Here we show the consistency between (|C4[) and (IC5I) . 
The starting point are the defining relations 

A = AJ Tl A — A Jt, (C6) 

where both A and A are in the Coulomb gauge. Whenever it is not confusing we use the compact notation 
A = AJt A 1 = Ay J 2 T . Moreover, the photon fields are related by the canonical transformation T such that 

A = TA, A f = A f T f (C7) 
The invariant object is basically the electromagnetic energy written in terms of the transverse current 

E=\J d 3 xd 3 y J^A^ix - y)4(y), (C8) 

where the sources are real. In momentum space this implies 

(4(*0)* = 4(-fc), (C9) 

and Eq. (|C8|) translates into 

£ = i J d 3 k J l T (~k)A i:j (k)J 3 T {k) = i y d 3 fc 4 A J T . (CIO) 

Reality of E further demands 

At(fc) = A(fc). (Cll) 

Expressing i? in terms of the fields yields 

E =^jd 3 k f T AJ T = ^Jd 3 k At (A" 1 ) A, (C12) 
in such a way that the equivalent description in terms of the barred quantities requires 

E=- I d 3 kA^ (A" 1 ) A. (C13) 



2 . 

Inserting the transformation (|C7[) we obtain the relation among the two propagators 

A = TATK (C14) 
The general structure of a Coulomb gauge propagator in our parity violating theory can be written as 

A = AiK + A 2 S, A^ A X K + A 2 S, (C15) 
where we have introduced the hermitian matrices 

K = [Kij] = [|k|% - kikj\ , S = [Sij] = [te"^k m ] , (C16) 

with the following properties 

S 3 = |k| 2 5, S 2 = K, (C17) 

in such a way that each propagator can be written in terms of the matrix S only. 
The canonical transformation T has the form 

T = al + [3S = T\ (C18) 
in momentum space, with real numerical coefficients 



^1 + ^(1+45 



> = = ^—^ , 1 a = -g 2 \k\ 2 . (C19) 
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Substituting in the relation (|C14[) and after some algebra we obtain the following conditions among the components 
of the respective propagators 

1 



A! 



A, = 



1- V|k| 
1 



[Ai + 2 5 A 2 ] , 
A 2 + 2 5 |k| 2 A 



From the expressions (|C4[) and (|C5|) together with the definition (|C15[) one can read 

k 2 a _ 2gk 2 



Ai 
Aj 



|k| 2 ((fc2)2_ V fc4|k| 2 ) : 

(fc 2 -V|k| 2 fcg) 
|k| 2 ((fc2) 2 _ W|k| 2V 



2^2 _ A n 2h4 |H 2n I ' 
2 



((fc 2 ) 2 -4 5 2 fc 4 |k 



A 2 



2.9|k| 



((fc 2 ) 2 ~4 5 2 fc 4 |k| 2 )' 



and verify that they satisfy the relations (|C20|) and (IC21I) . 



(C20) 
(C21) 

(C22) 
(C23) 
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In this Appendix we write the general form of the remaining LIV contributions to the electron self energy which 
final results are presented in Section VII, according to our general scheme of calculation 

(fc 2 -2 3 2 fc 4 |k| 2 /fc 2 ) 



A = 4ie 2 7 



d 4 fc 



(2tt) 4 [ fe 2 _ m 2] ^( fc 2)2 _ Ag 2\ k f k 4^ 



~ 2ie 2 
3 y 



c = —(gm) 



D + E = 32ie 2 m 



d 4 k 



ikrfc. 



2 1.2 



(2tt) 4 [ fe 2 _ m 2j /( fc 2)2 _ 4ff 2| k |2 fc 4 



ikrfe 



2 J.2 



d 4 fc 



(2tt) ( fc 2 _ m 2)2 ^ fc 2)2 _ 4g 2| k |2 fc 4 



J(k), 



d 4 k (fc 2 + |k| 2 /3) (fc 2 - 2g 2 fc 4 |k| 2 /fc 2 ) 
(2tt) 4 (fc 2 -m 2 ) 3 ( (fc 2) 2 _ V | k |2 fc4 j 



J(k), 



D = F 



E 



4ie 2 r d^k 
~ 9 J (2^F 



4ige 2 f d A k 



1 



4fc 2 



(2^) 4 



(fc 2 — m 2 ) 2 (fc 2 — m 2 ) 3 
3 ^ ||k| 2 



2 1.2 



((fc 2 ) 2 -4 5 



2 lkl 2 fc 4 



(fc 2 — to 2 ) 2 (fc 2 — to 2 ) 3 



ikrfc 



2 ^.2 



(fc 2 ) 2 -4 5 2 |k| 2 fc 4 



j{k). 
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